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Abstract 



This paper discusses the feasibihty of testing Bell's inequality with the charmonium decays rjc — > AA, 
Xco AA and 3/ip ^ AA. We develop a new formulation of the Bell's inequality represented with par- 
ticles' orientations which can be measured straightforwardly in experiment, and seek if it violates in these 
decay processes. It is illustrated that the rjc channel and the XcO channel maximally violate Bell's inequality, 
while the J/V' channel gives no inconsistency with it. Simulations dealing with experimental aspects are also 
implemented. The expected statistical fluctuations and achievable significance are estimated as a function 
of sample size. 



! 1 Introduction 

> 

With no doubt. Bell's inequality (BI) has played a significant role in modern physics. It establishes a 
clear discrimination between quantum mechanical and classical view of nature. Since when it was first 
I shown to be violated in a photons system by optical experiments T , what if with other type of particles or 
^ ■ systems have long been people's interest. Cases involving massive particles are in particular curiosity, since 
«) I generally a massive particle displays more classical properties. At present, some of the experiments have 
■ been implemented dealing with proton pairs or K^K'^ and B^B^ oscillations [3][1], however the variety of 
the experiments is still limited. 

In this context, charmonium decays cc AA have been anticipated to be one of the testing channel of BI 
[5] [6] in a highly exotic system: one involves high energy massive unstable fermions experiencing parity- 
rx] . violating weak interaction. In particular, rjc — > AA and XcO ~^ AA are exact realizations of Bohm's type 
^ I of EPR experiment [7^ in that these involve a spinless particle in the initial state decaying into two spin 
" " " one half fermions with opposite spins. In quantum mechanics (QM) this is known as an entangled state, 
leading to a strong correlation between the two fermions, that violates BI. J/-0 AA is similar, and is 
thought to be the most promising channel because of the abundant statistics. However this channel has 
the significant difference of having a spin 1 initial state. This allows involvement of relative orbital angular 
momentum between the generated hyperons A and A, which disturbs the entanglement and weakens the 
particle correlation. It has not been clear if this weakened correlation is nonetheless sufficient to violate BI 
(We discuss and give a conclusion to this in chapter 4). 

Practically, this test seems to include a complication that we have no means to directly measure the spin 
of decaying hyperons. For an indirect approach instead, it has been suggested by N. A. Tornqvist that the 
association between A (A) spin and the decay distribution of its weak decay A — ?> p7r~ (A -> p7r+) is useful 
for inferring the hyperon spins within the classical picture [5j . The distribution of this decay product proton 
(anti-proton) is known to be 

—— oc 1 -f aAfi ■ s — — cx 1 - a^n. • s' (1) 
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in the A (A) rest frame [5]. n [n') is the unit vector of the outgoing proton (anti-proton) direction and s (s') 
the polarization vector of the A (A) as shown in Fig. 2, which are ah defined in the A (A) rest frame, and 
aj^ are the decay parameters with experimental values of q;a = 0.642 ±0.013 |^ and ofA = if ignoring the 
small effect of CP violation. Due to the parity- violating nature of the weak interaction, the outgoing proton 
(anti-proton) prefers to fly along (against) the polarization of its parent A (A). In this sense, the hyperon's 
decay is its own polarimeter. It is worth pointing out this prescription is not valid if in QM in that it ignores 
the contribution from interference between the two spin states. In QM, the entire process cc — ?► AA — > pTrpTr 
should be treated as a coherent process in which the spins of the intermediate state (hyperons in this case) 
cannot in principle be well-defined. 

While Tornqvist suggested to test QM with the use of this tool [5] [6], S. P. Baranov extended it for testing 
local realistic theory (LRT) by reformulating BI with respective to the orientation of decay product p and 
p in the final state [T^] with the help of ([Ij. We basically follow this line but with attaching further and 
comprehensive discussion. To be specific, this paper consists of three sections: 



• The development of BI into a representation with experimentally measurable observables, using LRT- 
based arguments, where the processes cc — > AA and AA — )■ pTrpTr are considered separately. 

• The evaluation of whether cc — > AA — > pTrpTr has sensitivity of testing this newly derived BI by QM- 
based calculations. 

• An analysis of achievable significance of such a test and its experimental feasibility. 




Figure 1: An overview of the cc — > AA — > p7r^p7r+ process. First, a charm meson decays into a AA pair. 
These travel back-to-back in the meson's rest frame and decay into p7r^ and p7r+ respectively. We measure 
the orientation of p and p as the testing variables of BI. Experimentally, the charm mesons are supposed 
to be produced in e^e^ collisions. J/i/' is obtained directly while rjc and xco are generated via the decay of 
J/i/j and ?/;' respectively. 
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Figure 2: The relation of unit vectors defined in the A (A) rest frame and their projections, n, (n') is the 
orientation of outgoing proton (anti-proton), s (s') the polarization of A (A). Their projection onto a guide 
axis is labeled by the attached index. 



2 Momentum representation of Bell's inequality 

In realistic theories, a particle's spin is a definite physical quantity and normally treated as a 3 dimensional 
continuously valued vector just like classical angular momentum. Consider a two particles system with their 
spins of one half. The polarization vector s and s' follow the algebraic condition: 

\{SaS't)-{SaS',)\ < l + {Sbs',). (2) 

where a, b, and c are arbitrary unit vectors "guide axes" . s and s' have the norms of 1 and Sa is its 
projection onto a i.e. Sa — s ■ a; s^^ is that of particle 2 onto b; Sb, s'^ defined similarly (Fig. 2). Here 
s and s' are quoted as the realistic value of polarizations and considering the case in that they distribute 
probabilistically, and the ensemble average is weighted by the probability density function of them. These 
are seemingly the striking deferences from the conventional BI [TU] which has the exactly same form as ([2]) , 

I ( TOa "T-f, ) - ( "T-a ) I < 1 + ( TOf, m'^ ) (3) 

with nii and m[ being the result of measuring Si and s[ (i=a, b, c). ^ is valid no matter if (m^) is 
continuous or discrete providing — 1 < rrii {m[) < 1, or if the realistic values Si (s[) are different by each 
measurement, though the original Bell's discussion [10] takes discrete ones assuming measuring the same 
spin by a Stern-Gerlach type of experiment. Thus, since the both inequalities ([2|) and ([3]) has the same 
essence that they generally hold in (local) realistic view of physics and the constraints are due to the realistic 
interpretation of spin, we call the inequality ([2]) BI as well and use it to test LRTs. 

In extending to a relativistic case, in which the particles 1 and 2 belong to different respective frames, the 
CHSH (Clauser-Horne-Shimony-Holt) version of BI gives a more appropreate description 

\{sas',) + {sas'^) + {s,s[)-{s,s'^) \ < 2 (4) 

Here four guide axes are involved, two of which (a, c) are used for particle 1 and the other two (b, d) for 
particle 2. Each set of two guide axes is defined independently in their own frame. Recall that the original 
version of BI ([2]) has b in common to both frames. This leads to a confusion when in a relativistic case, as 
pointed out by |12j, while the CHSH inequality (j4]) remains well-defined with no such ambiguity. 
At this point, it is clear that (HJ is not capable of being tested by experiment because it is denoted by 
the realistic values instead measured values, that can yield different values due to the disturbance by the 
involvement of hidden variables, generally in LRTs. However we can translate it into a testable one, assuming 
the angular distribution ([T]) for the decay A — > p7r~ (A — > p7r+). If the two decays are independent, the 
correlations of hyperon spin and the proton (anti-proton) orientation can be tagged as 

( 4 ) = — ( (n • a) (n' . 5) ) ^— ( n', ) . (5) 
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Table 1: Kinematics of A and A from each charmonium decay, u is the fraction that the two decays are in 
space-like configuration: oj = 2 dti J^^^'^ dt2 -e^~ - = /3. 





Vc 


XcQ 




/3 


0.663 


0.757 


0.693 


7 


1.329 


1.621 


1.408 


OJ 


0.663 


0.757 


0.693 



n (n') is the traveling directions of the proton (anti-proton). The equation was originally provided by S. P. 
Baranov [12| . which however takes a different picture of classical spin reality from one in this paper. The 
derivation for our spin interpretation is given in the appendix, yielding the same result as |12j . Equation (|4]) 
can therefore be written as 

\{nan',) + {nan',} + {n,n',)-{nM\ < (6) 

We call this "the momentum representation" of BI. Here n and n' are measured values, in addition, have 
corresponding observables in QM. Therefore ([6|) is eventually an inequality which can be evaluated by ex- 
periment and by QM. 

Several comments should be added to this reformulation: 

• On deriving ([5]) , we assumed that the decay of the first hyperon only depends on the polarization of 
that particular hyperon. It seems reasonable according to the kinematical properties of A and A from 
each meson decay, as shown in Table 1. The two decays have a space-like separation for 66 % ^ 76 % of 
events, across which no interaction can act. Experimentally, these space-like events can be selectively 
extracted, which realizes a complete isolation of the decays. 

• In a conventional tests with direct spin measurement, guide axes play much important physical role 
in that only one direction of spin can be chose to be measured. In our test, however, guide axes are 
just arbitrary unit vectors as seen in ([2]). This cause no problem because momentum components can 
be simultaneously determined within one measurement, and even has an advantage in that we will not 
suffer from the "free will" problem in choosing guide axes. 

• On the other hand, it has no scheme comparable to " delayed-choice measurement" in optical experiment 
which prevents the hyperons from deciding the way of decaying in advance, or exchanging information 
each other just after their production so that they can have a larger correlation in their orientations. 
These are expected to be the loopholes of this test. 

• The development from (jl]) to ([6|) proceeds based wholly on the classical picture. Although ([T|) can be 
derived by a standard QM calculation, we accept ([T]) as just an experimental fact without assuming 
any background theories underlying it. 

• There has been a claim that BI with respect to commuting variables can not violate even in QM due 
to their incoherency [T3J. However this is not the case because the inequality ^ is essentially a BI 
with respect to spins. We just transform and represent it with n and n' , as a result, the upper limit 
of the BI is lower by factor of Hi^i -L) than that of a BI which naively takes variables as n (n'): 
\{n'an'ij) + {nan'^) + {ncn'^) — {ncn'^)\ < 2. Thus the BI ^ does violate, as we show in the section 4. 
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• It may seem to be against the intuition that © is testable inequahty regardless of the magnitude of 
QfA which represents the extent of the correlation between s and n. Theoretically this is true since the 
upper limit (the right-hand side of ^ drops accordingly. The conceptional essence of the test is not 
how much n {n') acts as the polarimeter of hyperons but that the correlations (sis'j) and {riin'j) 
have an exact relation ^ thus we can set the upper limit on ( rii n'^ ) as (jH])- Practically, however, the 
small correlation amplitude of ( ni n'^ ) is difficult to measure precisely, therefore we at last prefer the 
decay channel with as large aA as possible. 



3 Bilinear expression of the momentum represented BI ([6]) 

For later analysis, it is convenient to transform the BI ([6]) into one written with a "correlation matrix" C. 
C is a 3x3 real valued matrix given by the correlation amplitude, 

Cij :== {n,nj) {i, j ^ 1,2,3) 

where the indices i, j label x, y, z components of a vector in Cartesian coordinate, {uau'^) can be written 
in a simple form of bilinear using matrix C as 

3 3 

(riaO = Y^^a,hj{n,n'j) = C b . (7) 
i=i 

We now define Q as the left-hand side of which can be written as a sum of bilinears: 

Q^\a^C{b + d) + c^C{b-d)\. (8) 

The advantage of this representation can be seen in that the physical part {C) and physics-independent part 
(guide axes) are well separated. 

Next we try to specify the maximum value of Q and to get rid of guide axes out of the in equality since 
they give no physical importance in our formulation. It is easy to show that the maximum value of Q is 
same as that of aFC{b + d) + (F C{b — d), thus considering the case of Q = aFC{b + d) + (FC{b — d) 
is sufficient. The method of Lagrange multipliers (MLM) can be utilized. With the constraint conditions 
a^a = 1, b^b = 1, c = 1 and dF d = 1, a scalar function L can be constructed using four multipliers ^a, 
Cb, and . 

L = a^C (b + d) + c^C {b - d) 

- {a^a - 1) - {b^b - 1) - (c^c - 1) - \u {<fd - 1) 

Setting all the derivatives to zero, 

(8a) 
(86) 
(8c) 
(8d) 

(a — a, 6, c, d). (9) 

Mutiplying , c'"{b, d) from the left (right)-hand side of (8a) ^ (8d) respectively and using the constraint 
conditions (9), the multipliers can be written as 

- a'^C{b + d) = c'^C{b - d) 

ib = {a + cfCb ^d = {a- cfCd (10) 



dL 
doT 


= 


^ Cib + d)-^aa^ 


dL 
db 


= < 


=^ (a + cfC - ibb^ = 


dL 
d^ 


= < 


=^ C(b-d)-ecc = 


dL 
M 


= < 


=^ {a- cfC - = 


dL 
d^ 


= 4= 


=;> (constraint conditions 
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Substituting these back to (8b) (8d) gives 

^ _ C^{a + c) _ C^{a + c) 
^ (a + c)TCb ~' P 
,^C-ia^^,C^a-c^_ (11) 
{a-cYCd o 

Note that 6 and d are parahel to Cf^ (a + c) and [a — c). Normahzation factors p and cr are chosen to 
satisfy h^h = 1 and d^d = 1 i.e. 



p=^{a + c)C'^C{a + c) (j=\J{a~ c)C'^C{a - c) (12) 

Now we define a symmetric matrix S* as S* := C'^C = CC^ . With (|10p pT|) . we eliminate b and d in (8a) (8c), 
giving 

{a + p)Sa + {(J — p)Sc = /ia 
(ct — p)Sa + (cr + p)5'c = vc 

where 

^ := ((T + p)a^Sa + (ct — p)a^Sc 
V := (o- + p)c^Sa + (o- - p)c^Sc . 

These can be expressed as eigen-equations for the vectors a and c 

[-4/9(tS'^ + (cr + p)(^ + i^)5] a = ^i^a 
[^paS"^ + {(T + p){p + v)S\c = pz^c. (13) 

It can be easily verified that the matrices appearing in the left-hand sides of ([13]) ■izApaS'^ + {a+ p){p+v)S 
have identical eigen-vectors to those of S as 

[±ApaS'^ + {a + p){p + v)S] Vi = [±4,pa\j + {a + p){p + v, 

with Ai, Vi being the eigen- values and eigen-vectors satisfying Svi = XiVi. Since any 3-dimensional matrix 
can only have three independent eigen-vectors at most, Vi{i = 1,2,3) give a full description of all solutions 
of dUl). Therefore, 

a = ±Vi 

c = ±Vj 

\v,\^l (i,j = 1,2,3) (14) 



are required and we see all of these satisfy ([T3| . Note that the norms of a and c are confirmed to be 1. 
Substituting these into (ITTtp^ . 6, d and all other coefficients are determined. Setting (a, c) — {vi, Vj) for 
simplicity. 



{\) a = c{i= j) 



p = 2VA, (7 = 



d — (arbitrary unit vector) 



Q = 2a' Cb = 2VA, (15) 
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(ii) a^c{i^j) 



P = (T = ^/Xi+Xj 



d 



^ Xi + Xj -^Aj + Xj 

Q = 2 ^A, + Xj (16) 

In deriving we used the fact that S = C^C is a symmetric matrix with orthogonal eigen- vectors 
{vj ■ Vj = Sij). The same analysis can be apphed for the case of (a, c) = {vi, —Vj), {—Vi, Vj), {—Vi, —vj). 
The maximum Q value comes from (jl6p when we set A^ and A^ to the largest two eigen-values of S. One 
point to be noted is that solutions of MLM generally include local maxima where the real maximum value 
is given not at extrema but at the boundary of parameter space. However this is not the case here, since 
the parameter space of a, b, c, and d are independent "spheres" that have no boundary. 



With all these considerations, we finally reach the form. 



Qmax = 2yA7TA^ < (17) 

where Ai and A2 are the largest two eigen-values of C^C. For later convenience, we define C and Q by 
dividing by the right-hand side of (|17p 



C.,^{n,n'^)-^ (z,j = 1,2,3) 



Qmax = 2-\/ Ai + A2 (Ai, A2: the largest two eigen-values of C'^C) 

Qmax < 1. (18) 

This psp is our target expression. The classical limit of Qmax is Qcl = 1 while the quantum limit reaches 
V2. 



4 Quantum mechanical calculation of C and Qmax 

In this section, we perform a QM-based computation of the correlation matrix C and Qmax for each channel 
T]c — )■ AA — ^ p7r~p7r+, XcO — ^ AA — ^ p7r^p7r+ and J/V' — J' AA p7r~p7r+, to determine if the BI (fTS]) is held 
or violated in QM. 

The matrix elements for each channel are as below, according Feynman's rules and the effective Lagrangian 
prescription. 

Mrj, = 7V4a(Pa,sa,Pp,Sp)w(Pa:Sa)7^«(Pa>sa)-^a(Pa>sa>Pp.Si3) 

■^Xco = MA{PA^SA,Pp,Sp)u{pf„SA)v{pj^,Sx)Mx{pj^,Sj^,Pp,Sp) 
Mj/^ (X MA{pA,SA,Pp,Sp)u{pA,SA)e^ 

Ma{Pa,sa,Pp,Sp) = u(Pp,Sp) (1 + ca7^)'«(pa>sa) 
Ma{Pa, sa, Pp, Sp) = v{pA, sa) (1 - Ca7^) v{Pp, Sp) 

sa denotes the helicity of particle A, u(pj^,sa) and v{pa,sa) the 4-spinor and the polarization vector 
of the meson J/\['. Momenta appearing here Pa, Pa are all defined in the decaying meson rest frame. Ma 
and Ma are the matrix elements responsible for the A — > p7r~ and A — p7r+ sector which gives ca = ca 
providing CP conservation argument. The distribution of A decay and the decay parameter a a in ((T]) is 



^'(PA)Sa)>'a(PA>sa>Pp> 
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associated as 

(19) 

-2|p|cA . . 

with which the parameter ca is derived to c\ — C]^ — —6.79 ± 0.18. E, p and rup are the energy, momentum 
of outgoing proton in the A rest frame and its rest mass respectively. 

■Mj/^ includes an additional parameters associated with the form factor of J/tp- This has been ex- 
perimentally determined from the unpolarized A distribution of e+e^ — > J/^ — )■ AA where e+e~ act as 
chiral fermions leading J/^E* polarization parallel or anti-parallel to the beam axis according to the helicity 
conservation in a high energy system: 

2^ \Mj/^\ =1 + Fcos^0 



(l + r2)(A/2+4m2)-8Mmr ' 2^a^ + l 

M and m are the masses of 3/^ and A respectively. 9 is the angle between the e~^e~ beam axis and the direc- 
tion of A. The observed value of F is F = 0.65±0.11 [Hj which gives two possible solutions = 2.46±0.18 
and = 0.54±0.18. In our calculation here, both lead to identical results. (It probably depends only on F.) 

The angular distribution for the coherent process cc AA — !■ pirpi: is therefore obtained, assuming 
unpolarized initial and final states, 

2 

/ da 



\dflAdQpdQp 



SA.SA 



lA^ol'. (22) 



dflp {dflp) is the solid angle element in the proton (anti-proton) momentum space in A (A) rest frame. 
The average is taken over the spins of particles in initial and final states. Note that the sums over s\ and 
S]y are taken before squaring, which leads to the terms representing quantum interference between the two 
intermediating spin states. The correlation matrix C is calculated as 



Cij — ( fii n'j ) ^— 2" 



9 

2 
A 



[PtP'-i) 9 



\P\' 2ai 
9 



2«ilpP 



J dQAdfip dflp Pi p'j Prob(p, p' , p^) 

Jdn^dnpdnpp^p'^ (^^J^^^_/^to 



/ dQAdflpdQp Pi p'j \M 



|2 



01 



2"aIpP jdnAdQpdrip \Mo\ 



(23) 



p {p') is defined as the proton (anti-proton) momentum in the A (A) rest frame. (Recall that momenta Pp 

[pp) used in computing |A^o|^ above is one in the decaying meson rest frame.) The result of Cy and Qmax 
for each channel is shown in Table 2. The off-diagonal components are always 0, reflecting the symmetry in 
the processes, thus Qmax — + following (fT8l). 

Systematic uncertainties are estimated by shifting the parameters used in the calculation {e.g. particle 
masses, ca, etc.) within their Ict experimental uncertainties. The main contribution is from the uncer- 
tainty of while those from the other parameters have almost no effect on C or Qmax- One sees that Qmax 
in the rjc and the XcQ channel well surpass the classical limit Qcl = 1 and even reach the quantum limit \/2. 
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Table 2: Cij and Qmax for each channels are calculated as below. The uncertainty is dominated by those 
of measured parameters, especially a*. Values for r/c, XcO channels have only trivially small computation 
uncertainty since they are independent of measured parameters, in contrast to the i /ip channel. The classical 
limit for Qmax is Qcl = 1- 





Vc 


XcO 




Cll 


0.500 


-0.500 


-0.274 ±0.008 


C22 


-0.500 


0.500 


-0.177 ±0.024 


C33 


-0.500 


-0.500 


0.404 ±0.033 


Off-diagonal components 











Qmax 


1.414 


1.414 


0.976 ±0.046 


Qcl 


1 



Table 3: Some characteristics of the Qmax distribution 





Tic 


channel 






XcO 


channel 




n 


mean 


RMS 


significance 


n 


mean 


RMS si; 


jnificance 


100 


1.619 


0.653 


1.33 


100 


1.617 


0.652 


1.33 


300 


1.477 


0.403 


1.55 


300 


1.479 


0.405 


1.55 


500 


1.454 


0.316 


1.77 


500 


1.456 


0.314 


1.77 


1000 


1.432 


0.223 


2.20 


1000 


1.433 


0.225 


2.20 


2000 


1.423 


0.159 


2.86 


2000 


1.423 


0.159 


2.86 


3000 


1.420 


0.131 


3.43 


3000 


1.421 


0.131 


3.38 


5000 


1.417 


0.100 


4.24 


5000 


1.418 


0.101 


4.28 



In contrast, the J/'0 channel gives no significant excess in Qmax, therefore, is insensitive to test BI. The 7]c 
and the XcO channels conserve their entanglement throughout the process whereas the J/?/) channel loses the 
spin correlation in J/ip ^ A A where relative orbital angular momentum between the two hyperons dilute 
the spin correlation by some fraction. 



5 Estimation of necessary event number and the corresponding 
significance 

For the rjc and XcO channels, the BI is violated with the large Qmax values which however experience statistical 
fluctuations with a limited number of events. In this chapter, we show how the components of C and Qmax 
fluctuate statistically. Using MC simulations based on matrix (fT9|) . the distributions of Cn, C33 and Qmax 
are calculated as Figure 3, where we set the guide axes to the configuration that gives the maximum Q value 
as we discussed in section 3 i.e. Qmax — 2-\/ + Cfg. 

Cii and C33 fluctuate according to Gaussian distributions for all n. On the other hand, the Qmax distribution 
has a slight positive bias and can be well approximated by a Gaussian distribution providing n ^ 1000. The 
corresponding mean, RMS and significance are given in Table [3] The significance is calculated as the 
corresponding deviation in a Gaussian distribution of the p- value that Qmax fluctuates under the classical 
limit Qcl = 1- The two channels have the same n dependency and we see about 2000 events are sufficient 
to announce evidence of the BI violations. 
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Figure 3: Distributions of Cn, C33 and Qmax in the r/c and Xco channels. 
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Figure 4: Achievable significance with respect to number of events is estimated by MC simulation. We use 
a large number of samples where MC fluctuation is negligibly small. 
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Table 4: Branching fractions of each decay in J/i/j — > 777c ; ?7c ^ AA — )• pTrpTr and 

V'' 7Xco; XcO ^ AA ^ pTrpTT. 



channel 


branching fraction 


J/V' ^ + 7 
■0' ^ XcO + 7 
r]c^ A + A 
XcO ^ A + A 

A ^ p + TT" 

A-> p + 7r+ 


(1.7 ±0.4) X 10-2 
(9.7 ±0.3) X 10-2 
(1.41 ±0.17) X 10-3 
(3.3 ±0.4) X lO-'' 
(6.39 ±0.05) X 10-1 
(6.39 ±0.05) X 10-1 


V'' ^ 


7'?c; r?c AA ^- pTrpTT 
7Xco; XcO ^ AA ^ PTrpTT 


(9.8 ±2.6) X 10-s 
(1.31 ±0.16) X 10-5 



The number of experimentally available events and the measurement feasibility are also studied. The 
branching fractions of each relevant decay are listed in Table 4, assuming that rjc and XcO are all produced 
via J/V' — >■ f?c + 7 and tp' — >■ XcO + 7 respectively. In the BES3 experiment, 1 x 10^ of J/ip and 4 x 10^ of 
■0' are planned to be produced by the end of 2012 [15]. Using only the space-likely separated events which 
account for around 70% of all events (Table. 1), the event yields are about 6500 x e and 4000 x e for the rjc 
and the XcO channel respectively, with the event acquisition efSciency e. The r]c channel has enough yield 
even with conservative selection of events (efficiency e ~ 0.3) while the XcO channel is available provided a 
looser selection with e ^ 0.6. 

6 Conclusion 

Charmonium decays rjc — > AA, XcO — > AA and J/'0 — > AA are possible probes of Bell's inequality, due to 
the strongly correlated spins of the hyperon pair. As the hyperons undergo decays in which the angular 
distribution is associated with the hyperon spins. Bell's inequality can be developed into a new expression 
in terms of a momentum correlation of the decay products p and p. By QM calculation we find that the 
correlation in the J/^ channel is not strong enough to test Bell's inequality, while it is sufficient in the other 
two channels. The violation can be experimentally confirmed with around 2000 events for each channel; this 
number of events would have already been produced at BES, and this experiment is therefore now feasible. 
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Appendix: The derivation of equation ([5]) 

We start with the angle distribution of the protons (anti-protons) from A (A) decays (Same equation as ((T|) 
in the section 1): 

P{n\s) ^ 1 + UATi ■ s P{n'\s') = 1- aj^n' ■ s'. (24) 

n(n') is an unit vector of the the proton (anti-proton) orientation. P{n\s) {P{n'\s')) indicates the condi- 
tional probability density of n {n') with given polarization of the hyperon s (s'). Note that these distribu- 
tions are normalized to 1 with n (n') integration over solid angle / (ir2„/47r (J diln' / 4:Tt) . The correlation 
amplitude ((n • a)(n' • b)) can be calculated in terms of these distribution. 

((n.a)(n'.6))= J ^^?^(n • a)(n' • 5) P(n, n') 

-^^^^{ri ■ a){n ■ b) P{n, n\s,s) P{s, s ) 

P{x) represents the probability density function which a set of variables x follow. When the two hy- 
peron decays are isolated in space-time, which we confirmed in the chapter 2, the joint probability density 
P{n\s)P{n'\s') should be identical to the combined probability density P{n,n'\s, s') since they have no 
correlation in between, according to the instruction of the locality principle. Using , therefore, 

{{n-a){n'-b)) 

= I ^^^^iri-a){n'-b)P{n\s)P{n'\s')P{s,s') 

= J ^^^^(-•«)(-'-b)(l + «A--)(l-«An'-.V(.,A 

dflndQn' integration can be performed using spherical polar coordinates {0 ,(/)), {9' , (j)' ) with s,s' being the 
zenithes: 

{{n-a){n'-b)) 

/ d cos 9 d(t) d cos 9' d(j)' drts d^ls' , . „ ■ n ■ i 
(ai sm 9 cos m + a2 sm 9 sm m + cos 9) 
An 47r 47r 47r ^ v o ; 

X {bi sin 9' cos0' -I- 62 sin 6*' sine/)' + 63 cos6'')(l -I- a\ cos6')(l — cos 9')P{s, s') 

Here cos(p and sin0 terms vanish in the integral, thus only cos 9 terms remain. As 03 = a • s = Sq and 
bs — b ■ s' — s'fj, we obtained the desired form: 

{{n-a){n'-b)) 



f dcos9 dcos9' ,^ ^, ^ ^, 

/ — — (1 -I- q;a cos t^)(l — cos ) cos cos 
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